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X-Ray View of the Shock Front in the Merging Cluster Abell 3376 with Suzaku* 
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Abstract 

We report on a Suzaku measurement of the shock feature associated with the western radio rehc in the merging 
cluster A 3376. The temperature profile is characterized by an almost flat radial shape with kT ~ 4 keV within 
0.5r20{) and a rise by about 1 ke V inside the radio relic. Across the relic region (0.6 - 0.8r2oo), the temperature shows 
a remarkable drop from about 4.7 keV to 1.3 keV. This is a clear evidence that the radio relic really corresponds to 
a shock front possibly caused by a past major merger. The observed sharp changes of the temperature and electron 
density indicate the Mach number M~ 3. The radial entropy profile is flatter than the prediction (r' ') of numerical 
simulations within 0.5r2oo, and becomes steeper around the relic region. These observed features and time-scale 
estimation consistently imply that the ICM around the radio relic has experienced a merger shock and is in the 
middle of the process of dynamical and thermal relaxation. 

Key words: galaxies: clusters: individual (Abell 3376) — galaxies: intergalactic medium — shock waves — 
X-rays: galaxies: clusters 
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1. Introduction 

Clusters of galaxies are believed to grow through gas ac- 
cretion from large-scale filaments and mergers of subclusters. 
Cluster mergers are the most energetic events in the Universe 
after the Big Bang with the total kinetic energy of the colliding 
subclusters reaching 10^^ ergs. The kinetic energy is converted 
to thermal energy by driving shocks and turbulence. Existence 
of mega-parsec scale radio emission (halos and relics) in merg- 
ing clusters indicates that those shocks and turbulence operate 
as main mechanisms of particle acceleration (see Ferrari et al. 
2008 and reference therein). Radio relics are considered to be 
the synchrotron emission generated through the interaction be- 
tween relativistic electrons, accelerated by a shock, and ampli- 
fied magnetic field in the intracluster medium (ICM). Although 
these radio features provide direct evidence of particle accel- 
eration, the detection of non-thermal (hard) X-ray emission 
caused by relativistic particles is still controversial (Nevalainen 
et al. 2004; Ajello et al. 2009). CuiTently, hard X-ray observa- 
tions set lower limits of magnetic fields in ICM (Wik et al. 
2009; Sugawara et al. 2009; Clarke & Ensslin 2006). Based on 
recent studies, it is believed that magnetic field of 0.1-10 /vG 
exists in ICM. 

Theoretical studies of cluster merger shocks provide rich in- 
formation about the magnetic field, turbulence, particle accel- 
eration and their time evolution (see e.g. Briiggen et al. 2011 
for a review). Merger shocks are characterized by a Mach 
number Al < 3, causing amplification of magnetic fields to 
Z? ~ 10 fiG and production of relativistic particles (Takizawa 
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& Naito 2000; Takizawa 2008). Also, turbulence is considered 
to give a significant pressure support in the outer regions, with 
around 10% of the thermal pressure at least for a few Gyr after 
a major merger (Parrish et al. 2011). Recent numerical simu- 
lations include interaction of merger shocks with the filamen- 
tary cosmic web structures outside of clusters, and the results 
indicate that shock fronts are more enhanced in the filament 
direction (Paul et al. 2011). 

Though past X-ray observations have shown images and 
temperature structures of cluster mergers, there are only few 
clusters for which clear evidence of the shocks has been 
obtained (1E0657-56, A520, and A2146: Markevitch et al. 
2005; Clowe et al. 2006; Russell et al. 2010). The low X-ray 
surface brightness in the outer regions has hampered precise 
measurements of gas temperature and density in the pre-shock 
region in particular. Recently, existence of shock fronts with a 
Mach number of about 2 was confirmed in the radio relic region 
of A 3667 (Finoguenov et al. 2010; Akamatsu et al. 2011b). 
Additionally, a remarkable radio relic indicated a shock front 
with a Mach number of 4.6 based on the measurement of ra- 
dio spectral index (CIZA2242: van Weeren et al. 2010). Thus, 
radio relics are the good probe to identify merger shocks and 
combination with X-ray data will yield important information 
about the gas dynamics associated with cluster mergers. Since 
the radio relics are mainly found in the outskirts of clusters 
where the gas density is low, sensitive X-ray observations are 
needed, especially in the upstream side of the shocks. The X- 
ray data will allow us to estimate the parameters of the gas and 
the magnetic field, and enable us to look into the actual process 
occurring during the cluster evolution. 

In this paper, we present a new X-ray evidence of the 
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Fig. 1. X-ray image of A 3376 in tlie energy band 0.5-8.0 keV, after 
subtraction of the NXB with no vignetting correction and after smooth- 
ing by a 2-dimensional Gaussian with o" = 16 pixel =17". The large 
blue circle shows the virial radius of A 3376 with r = 34'. 6, and white 
dotted circles show the annular regions used for the spectral analy- 
sis. Small magenta circles show point sources which are detected by 
wavdetect (see text). The VLA 1.4 GHz radio image is shown with 
green contours. 

shock associated with the radio rehc in Abell 3376, a well- 
known merging cluster with irregular morphology at z = 0.046. 
Previous studies showed the global mean temperature to be 4.0 
keV and the existence of a pair of Mpc-scale radio relics re- 
vealed by 1.4 GHz VLA NVSS observations (Bagchi et al. 
2006). Another feature of A 3376 is the report of the hard 
X-ray signal with BeppoSAX/PDS (Nevalainen et al. 2004). 
However, Suzaku HXD observation gave an upper limit which 
did not exclude the BeppoSAX flux (Kawano et al. 2009). We 
need to await further sensitive observations in the hard X-ray 
band. 

We use Ho = 70 km s"' Mpc"', Qm = 0.27 and Qa = 
0.73, respectively, which give 54 kpc per arcminute at 
z = 0.046. The virial radius is approximated by r2oo = 
2.77h-^({T}/10keyy'^/Eiz) Mpc, where E{z) = (Qm(1 +zf + 
1 - Qm)'^^ (Henry et al. 2009). For our cosmology and red- 
shift, r2oo is 1.86 Mpc (= 34.'6) with kT = 4.7 keV. We em- 
ploy solar abundance defined by Anders & Grevesse (1989) 
and Galactic absorption with A^h - 5.8 x 10^*'cm"^ (Dickey & 
Lockman 1990). Unless otherwise stated, the errors correspond 
to 90% confidence for a single parameter 

2. Observations and Data Reduction 

As shown in Fig. 1, Suzaku carried out two pointing obser- 
vations of Abell 3376 along the merger axis in October and 
November 2005, designated as A 3376c, A 3376w for the cen- 
ter and the west relic. The observation log is summarized in 
Table 1 . All the observations were performed with either nor- 
mal 5 X 5 or 3 X 3 clocking mode. The combined observed field 
extends to the virial radius of A 3376 (34'.6 ~ 1.86 Mpc) indi- 
cated with blue circle in Fig. 1. 

The XIS instrument consists of 4 CCD chips: one back- 
illuminated (Bl: XlSl) and three front-illuminated (Fl: XISO, 
X1S2, X1S3) ones. Even A 3376 observations were car- 
ried out only 3^ months after the launch, the IR/UV block- 
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Fig. 2. The spectrum of the Q0551 field used for the background esti- 
mation, after subtraction of NXB and the source. The XIS BI (black) 
and FI (red) spectra are fitted with CXB + Galactic components (LHB, 
MWH) (apec+wabs(apec+powerIaw)). The CXB spectrum is shown 
with a black curve, and the LHB and MWH components are indicated 
by green and blue curves, respectively. 

ing filters had some contamination accumulated at the time 
of the observations. We included its effect and uncertainty 
on the soft X-ray effective area in our analysis. We used 
HEAsoft version 6.11 and CALDB 2011-06-30 for all the 
Suzaku data analysis presented here. We performed event 
screening with the cosmic -ray cut-off rigidity (COR) > 8 GV 
to increase the signal to noise ratio. We extracted pulse- 
height spectra in 10 annular regions whose boundary radii were 
2', 4', 6', 9', 12', 15', 18', 21', 24', 27' and 31', with the cen- 
ter at (06'>01'"00% -39°57'07"). We analyzed the specti-a in the 
0.5-10 keV range for the Fl detectors and 0.5-8 keV for the BI 
detector. The energy range of 1.7-1.9 keV was ignored in the 
XIS spectral fitting, because the response matrix around Si-K 
edge had some residual uncertainties' . In all annuli, positions 
of the calibration sources were masked out using the calmask 
calibration database (CALDB) file. 

3. Spectral Fitting 

3.1. Method of Spectral Fits 

In the analysis of A 3376 data, we followed our previous 
analysis for the Suzaku observation of A 3667 (Akamatsu et 
al. 201 lb), which shared many similar features but with much 
brighter surface brightness than the present object. The ba- 
sic method of the spectral fit is described in subsection 4. 1 of 
Akamatsu et al. (201 lb). The observed spectrum was assumed 
to consist of thin thermal plasma emission from the ICM, 
the local hot bubble (LHB) and the milky way halo (MWH) 
as the Galactic foreground components, cosmic X-ray back- 
ground (CXB) and non-X-ray background (NXB). The NXB 
component was estimated from the dark Earth database by the 
xisnxbgen FTOOLS (Tawa et al. 2008) and was subtracted from 
the data before the spectral fit. To adjust for the long-term vari- 
ation of the XIS background due to radiation damage, we ac- 
cumulated the NXB data for the period between 50 days before 
till 150 days after the observation of A 3376. 

' http://heasarc.nasa.gov/docs/suzaku/analysis/sical.html 
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Table 1. Suzaku observations of A 3376 and the background region. 



Target (Obs. ID) 


Position 


Start time 


Exposure* 


Exposure^ 




(R.A., Decl.) 


(UT) 


(ks) 


(ks) 


A 3376c (100034010) 


(90.56,-39.94) 


2005/10/06 14:46:08 


110.4 


78.0 


A3376W (800011010) 


(90.05,-39.98) 


2005/11/07 14:15:05 


120.0 


99.9 


Q 0551-3637 (703036020) 


(88.19,-36.63) 


2008/05/14 13:25:55 


18.8 


15.3 



*: no selection with COR2 t: COR2 > 8 GV 



The spectrum of the Galactic background emission can be 
represented by two thin-thermal plasma models, by apec with 
1 solar abundance. The estimation of the Galactic compo- 
nent and CXB is described in subsection 3.2. Because the 
Galactic and CXB components have almost uniform spatial 
distribution in the XIS field of view, we generate uniform aux- 
iliary response files (ARE) for those emission. As for the ICM 
emission, we generated the ARE assuming the jS-model surface 
brightness with p = 0.40, and = 2.'03 (Cavagnolo et al. 2009) 
using xissimarfgen (Ishisaki et al. 2007). 

Eor the spectral fits, we used XSPEC verl2.7.0. We car- 
ried out spectral fit to the pulse-height spectrum in each annu- 
lar region separately. In the simultaneous fit of the BI and FI 
data, only the normaUzations were allowed to be different be- 
tween them, although we found that the derived normalizations 
were quite consistent between the two within 15%. In the cen- 
tral regions within 25 .'2, free parameters were temperature kT, 
normalization norm and metal abundance Z of the ICM com- 
ponent. In the outer regions (> 25 .'2), we fixed the ICM metal 
abundance to 0.2 as reported by Kawano et al. (2009). 

The 0.1-2.4 keV luminosity of A 3376 {Lx = 2.5 x 
10'*'* erg s"*) is less than 1/3 of the A 3667 value {Lx = 
8.8 X 10"*^ erg s"*), based on the ROSAT X-ray-brightest Abell- 
type clusters survey (Ebehng et al. 1996). To examine the ef- 
fect of flux contamination from difl'erent spatial regions, we 
estimated mutual photon contributions among all the annuli 
(see Appendix. 1). The resultant photon fractions show that 
most of the detected photons come from the "on-source" re- 
gion and about 20% is from the adjacent annular regions. 
These flux contamination features are very similar to that in 
A 3667 (Akamatsu et al. 2011b), in which the effect of stray 
light is always smaller than the systematic error of NXB. 

Recent study of the stray light properties of Suzaku XRT 
(Takei et al. 201 1) showed that its effect was significantly lower 
when the observing roll angle was inchned by 45 degrees just 
as in the case of the A 3376 west relic pointing. Based on this 
estimation, we can expect that the effect of the stray light does 
not influence the result significantly. We will not deal with the 
stray light any further in this paper. 

3.2. Estimation of the Background Components 

To exclude the point source contamination, we use wavde- 
tect tool in CIAO package version 4.3.0^ to identify point-like 
sources. We adopt the value of the sigthresh parameter as 
5 X 10~* in wavdetect, which is the significance threshold for 
the source detection. As shown by magenta circles in Fig. 1, 
we detected 12 point-Uke sources and masked out the regions 

^ iittp://cxc.iiarvard.edu/ciao/ 



with 1 arcmin radius from the sources. The lowest detected 2- 
10 keV flux of these sources was 5c = 5 x 10"*"* erg cm"^ s"* . 
Since this flux limit is lower than the level of 5^ = 2 x 10"*^ 
erg cm"^ s"' by Kushino et al. (2002), we adopted the CXB 
intensity as 5.97 x 10"^ erg cm"-^ s"* sr * after these point 
sources were subtracted. We also estimated the fluctuation of 
the CXB intensity based on Akamatsu et al. (2011a). Table 3 
shows the resultant CXB fluctuation assuming the flux limit as 
5c = 5xl0-*'*ergcm-2s-*. 

Since A 3376 is located far from the Galactic center (/ = 
246.53°,/? = -26.29°), we did not consider contribution from 
the emission associated with the Galactic plane. To estimate 
the Galactic background, we use another Suzaku observation 
data for Q 0551-3637 (Observation ID: 703036020, Hereafter 
Q0551)at3°.8fromA3376 as shown in Table 1. We excluded 
the region of the main source Q055 1 with a radius of 3'. Using 
the data in the Q0551 field, we evaluated the model of the 
Galactic background described in Sec 3.1 with A^h = 3.2 x 10^** 
cm"^. Eor the fitting, we fixed the temperature of the LHB 
component to 0.08 keV. The spectrum was well fitted with the 
above model (x1 - 1-08 for 135 degrees of freedom). The re- 
sultant parameters are shown in Table 2. The temperature of 
the MWH is 0.27^qq3 keV, and the temperature and intensity 
are consistent with the typical Galactic emission. For the esti- 
mation of the systematic error of the background spectrum, we 
adopted the CXB fluctuation as 15%, and the OBF contamina- 
tion to be ±10%, respectively. 

In the next section, we also examined the effect of systematic 
errors on our spectral parameters. We considered the system- 
atic error for the NXB intensity to be ±4.5 % (Tawa et al. 2008), 
the fluctuation of the CXB as shown in Table 3, and the OBF 
contamination to be ±10%, respectively. The resultant param- 
eters after taking into account the systematic errors are shown 
in Table 2. 

3.3. Single-Temperature Fit 

The pulse-height spectra and the best-fit models for all the 
annular regions are shown in Fig. 3. The parameters and the 
resultant values are listed in Table 4. Considering both the 
statistical and systematic errors, the ICM signal was signifi- 
cantly detected out to 31' (1.67 Mpc ~ 0.9r2oo) for the first 
time. 

Fig. 4 shows radial profiles of ICM parameters (temperature, 
surface brightness, deprojected electron density, and thermal 
pressure). The temperature within 20' (1.2 Mpc corresponding 
to 0.6r2oo) of the cluster center shows an almost flat or rising 
profile with kT ~ 4.0 keV, which is very much different from 
those seen in the relaxed clusters. The increase of temperature 
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NOMINAL 


CXBMAX 


CXBMIN CONTAM1+ 10% 


CONTAMl-10% 


LHB 












kT(keV) 


0.08 (fix) 


0.08 (fix) 


0.08 (fix) 


0.08 (fix) 


0.08 (fix) 


norm* (xlO^^) 


11 3+^'* 


10.6!^:? 


11.4!3-2 


Q 9+3.3 
^•"^-4.5 


"•^-3.4 


MWH 












;tr(keV) 


rv 97+0.06 
^■^'-0.03 


n TO +0.06 
*'-^^-().()4 


0-28!°:°^ 


ri 97+0.06 
"■^'-0.04 


n 77+0.06 
"■^'-0.04 


norm*(x 10""*) 




4 6+'-'' 

^•"-1.4 


^•'-1.4 


r 9+2.2 






143/ 135 


147 / 135 


152 / 135 


145/ 135 


143/ 135 



*: Normalization of the apec component scaled with a factor l/400;r. 

^orm=-^^ ^nenndV I (An{\ +z^)Z)^) x 10"'^ cm"^ arcmin"^, where Da is the angular diameter distance to the source. 



Table 3. Estimation of the CXB fluctuation. 



Region 


0'-2.'0 


2.'0-4;0 


4r0-6.'0 


6:0-9:0 


12:0-15:0 


15:0-18' 


18'-21' 


21'-24' 


24: -27' 


27'-31' 


SRR' 


4.01 


3.35 


3.29 


3.40 


0.36 


0.76 


1.08 


0.95 


0.83 


0.66 


e,Suzaku 


12.4 


18.0 


30.3 


50.2 


12.6 


34.6 


60.7 


64.5 


67.0 


63.7 


o"//cxb1; 


47.1 


39.1 


30.1 


23.1 


28.2 


21.3 


20.6 


20.3 


20.8 


20.2 



* ■.SOURCEJiATIOJiEGi%] (See Akamatsu et al. 2011a Sec 3.4) 
t: [arcmin^] 

%: Sc = 5x 10"''' erg cm"^ s"' is assumed for all regions. 



Tables. The best fit parameters of 2 kT model with statistical 
and systematic errors (CXB fluctuations and OBF contamination) 

2kT 

kTugh 4.05+1.62+1.37 + 0.11 
kTiow 1.09+0.22 ± 0.25 ± 0.07 
Z 0.2 (fix) 

norm'ig^ 1 .56 + 0.6 ± 0.7 ± 0.3 
norml^ 0.12+0.03 + 0.8 ± 0.2 
X^/ d.o!f 329/306 

*: Same as Table 4. 

t: 10-40 keV power-law flux in units of 
lO"''' erg cm~^ s"' arcmin"^ 

just inside of the radio relic (18'-21', and 21'-24') compared 
with the level in the inner region is significant by about 4cr. 
The temperature of A 3376 suddenly drops at the radio relic 
region from 4.7 keV to 1.3 keV. These features strongly sug- 
gest that there is a shock heating taking place in the relic region. 
The surface brightness profile also shows more than an order of 
magnitude drop at the same region. 

Based on the same Suzaku data, Kawano et al. (2009) 
showed the temperature in the radio rehc region to be /cr a; 3.8 
keV. The region roughly corresponds to a combination of 
21' - 24' and 24' - 27' in our analysis, and our temperatures 
of 4.7 keV and 3.1 keV for the two regions indicate that our 
results are consistent with the previous study. 

We also calculated the thermal pressure profile, with the 
pressure defined as P = riekT. The deprojected electron density 
He was derived from the apec normalization as in Akamatsu et 
al. (201 la). As seen in Fig. 4, the deprojected electron density 
and pressure profiles show a remarkable drop across the radio 
rehc, which strongly suggests existence of the shock front. 



3.4. Spectral Features in the Radio Relic Region 

The radio relic region was mostly covered by the deep 
Suzaku observation (~ 100 ks). In the radio relic of A 3667, 
Akamatsu et al. (2011b) examined the possibility of multi- 
temperature plasma using 2kT model and several DEM 
(Differential Emission Measure: Kaastra et al. 1996) models. 
They confirmed that the emission measure distribution is char- 
acterized by two peaks around the radio relic region. The sharp 
drop of temperature in the radio rehc of A 3376 suggests that 
a multi-temperature model is needed to describe the spectrum. 
The data will allow us to examine the relative contributions 
from different temperature components. 

To examine the multi-temperature feature in the radio relic 
region (24'-27'), we included additional thermal component 
(2kT) to the single temperature model. We then obtained an 
almost acceptable fit with;^;^ d.o.f. = 329/306 (Table 5) com- 
pared with the single temperature case of 345/308. The rather 
low statics of the spectrum hampered us to distinguish be- 
tween the 2kT and IkT models. The hot component temper- 
ature was derived as kT^ig^ = 4.05 ± 1.62 keV, and the cool 
one as kTiow = 1.09 + 0.22 keV, while the common abundance 
was fixed to 0.2 solar. These values seem reasonable con- 
sidering the result of the single temperature fit, which gave 
kT = 3.11 ± 1.10 keV. These two temperature values are close 
to those observed inside and outside of the radio relic regions. 

The actual three dimensional structure of the shock front 
would be a "bow" shape. Considering that the observed data 
should contain both the upstream and downstream components 
with respect to the shock front, it is natural to expect that the 
data in the radio relic region of A 3376 is a mixture of high 
and low temperature ICM by the projection effect. The ratio 
of emission measures can be used to constrain the shock front 
structure, however we leave such a study for future works. 
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(a) 0'-2' 




Energy (keV) Energy (keV) Energy (keV) 



Fig. 3. NXB subtracted spectra in each annular region. The XIS BI (black) and FI (red) spectra are fitted with the ICM model (wahs + apec), along with 
the sum of the CXB and the Galactic emission (apec + wahs(apec + powerlaw)). The CXB component is shown with a black curve, and the LHB and 
MWH emissions are indicated by green and blue curves, respectively. The total background components are shown by the orange curve. 



6 



H. Akamatsu et al. 
Table 4. Best-fit parameters of the ICM 
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0'-2' 


2' -4' 


4' -6' 


6' -9' 


12' - 15' 


15' -18' 


18' -21' 


21' -24' 


24' -27' 


27' -31' 


NOMIMAL 


kT (keV) 
Z 
norm' 

0.4-10 keV 

Ar7d.o.f 


^■"^-0.08 

33+"-''3 

"■-^-^-0.03 

164.8!3| 
70 « 
322/307 


4 16+""^ 

n 9O+0.03 

108.3+;| 
46.28t;3 
314/307 


4.14+""' 

26+^-^^ 

"■^"-0.03 

75.0+J-J 
388/307 


4.17+""' 
25+^'^^ 

" -^-0.03 

64.7+JI 
26.25!i-« 
342/307 


4.08!!J:^? 

46.4+11 
19 04+" ** 
362/307 


4 17+"-^' 
n 1 s+0.06 

25 5+" ' 
10 92+"-28 
352/307 


4 3J+0.29 

^■°^-0.28 

11 +"06 
"■^'-0.06 

15 3+"-5 

•'•■^-0.5 
6 14+0.22 
"•^ -0.22 

348 / 307 


4.68!»:^« 

0.2 (fix) 

6 7+"' 
"•'-0.3 
9 79+0.15 

-0.15 

351/308 


3 11+11" 

■'■^^-0.78 

0.2 (fix) 
2 l+"-2 
77+"04 

"•"-0.04 

345/308 


1 34+0.69 

-0.39 
0.2 (fix) 

3+"-2 

13+""i 

-^-0.01 

361 / 308 


CXB MAX+NXB 4.5% RED 


kT (keV) 
Z 
norm* 

0.4-10 keV 

P/d.o.f 


4 1 5+0.118 
'^■•'-'-0.03 

164.8+^^ 
70.95!l:^8 

327/307 


n 9O+0.03 

108.3+1' 
46.89!°;ii 
317/307 


4 22+"-"'' 
0:26+11 

76.1!!:! 
'19 09+0.14 

■-'^■■J^-0.14 

391 / 307 


4 26+0.11 

^■^°-0.09 

65.4!li 
26.57!i« 
347 / 307 


^■^^-0.28 

o-2o!S:!!^ 

20 35 +"-55 
362/307 


4 56+"-^= 

^■^"-0.25 

19+0.06 

^^■^^-0.06 

26 2+" ' 
11 36+"-2* 
353/307 


5 34+0.36 

-0.28 

0.12!iJ;«6 
15.8!":^ 

6 45+"-2t 

"•^-'-0.24 

350 / 307 


5 57+"-" 

■^•■"-0.48 
0.2 (fix) 

7 n+0.3 

'■^^-0.3 

3 02+" " 

-0.17 

356/308 


5 35+'-''= 
■^■-^■^-1.00 

0.2 (fix) 

2 6+"-2 
^■°-0.2 

1 07+""* 

^■"'-0.06 

352/308 


9 97+1.65 
Z.Z/ j 

0.2 (fix) 

5+" 2 
'^■■'-0.2 

21+" "2 
-0.02 

373 / 308 


CXB MIN+NXB 4.5% ADD 


kT (keV) 
Z 
norm* 

0.4-10 keV 

Ar7d.o.f 


4.04+»-^ 

-0.03 

161.6!f;^ 
69.04!|-^^ 
316/307 


4 09+"u« 

-Q.Q9 

n 9O+0.03 

106.4+11 
45.49!«;|1 
314/307 


4 06+"-"' 

^■""-0.09 

26+"-"' 

"■^"-0.03 

72.8![-| 

31 25+"" 
-'^■^-'-0.11 

391 / 307 


4 09+0.09 

0.25!»|| 
63.4![i 
25.76!]-3^ 
339 / 307 


09+0.29 

22+"-"' 
''■^^-0.09 

42 7+2-3 
17 53+" " 
367 / 307 


3 99+U.21 

"■l»-0.06 
24 2+" ' 

350/307 


4 52+0.31 
^■■^^-0.28 

12+" "6 

"■'^-0.06 

14 0+"-^ 

^■™-0.19 

357 / 307 


4 00+"'" 

^■""-0.40 

0.2 (fix) 

5 7+"-' 
-'■'-0.3 

2 2Q+"-'2 
^■^^-0.12 

367 / 308 


1 55+"-26 

i.jj 21 

0.2 (fix) 

1 l+"-2 

41+0-02 

"■^'-0.02 

351/308 


0.87+;;- 

0.2 (fix) 
l+"-2 

04+" "' 
"■'-'^-0.01 

414/308 


CONTAMI 10% ADD 


kT (keV) 
Z 
norm* 

0.4-10 keV 

P/d.o.f 


^■''■'-o.os 

"■^ -0.03 

164.8+^;^ 

70 12+'® 
'"■^ -1.69 

323 / 307 


4 j2+<"W 
(1 90+8:8? 

108.3+]^ 
46 10+""*^ 
316/307 


4.12+"-"' 
26+^-^^ 

■^■^"-0.03 

75.0+[-| 
31 69+"^' 
388/307 


4 16+""' 
^■^ -0.09 

26+" "' 

^^■^"-0.03 

64.111;^ 
25 92+'-''' 
338 / 307 


3 96+"-2' 
21+" "' 

"■^^-0.08 

45.9+^-^ 
34 94+1.59 
^^■^^-1.59 
364 / 307 


4 20+"-'^^ 
^■^"-0.20 
n 1 s+0.06 

"■l**-0„06 

25 1+"' 

10 77+" 21 
lu. / /_Q 2i 

348/307 


A SK+0.29 

12+""'' 
^'■'^-0.06 
14 9+0.5 

6 03+" 25 
354/307 


4 73+0.48 
•'■^-0.46 

0.2 (fix) 

^;+o.3 
6.6_o3 

2 64+"'* 

-0.16 

363 / 308 


3 16+1" 

-0.80 

0.2 (fix) 

1 9+0.2 

A 79+0.04 
-0.04 

360/308 


1 37+1-''' 

i -J' -0.40 

0.2 (fix) 

2+"i 
'^■^-o.i 
n n^;+o.oi 
0.06_poi 

402/308 


CONTAMI 10% RED 


fcr (keV) 

Z 

norm* 

5 + 
0.4- 10 keV 

^-/d.o.f 


4 12+»i)8 


"■■^■-^-0.03 

161.6+fl 

69.34![-^5 
317/307 


4 18+""' 
-Q.Q9 

n 90+0.03 
"■^°-0.03 

106.4+jl 

45 65+" '*'* 
315/307 


4 1 5+0.1)9 
n 9zr+0.03 

73.9+;-[ 

31 53+"''" 

Jl.JJ 03^ 
390 / 307 


4 16+"-"' 

^■'°-0.09 

25+"-"' 

"■^-'-0.03 

64.7+1] 
25.92!i;^^ 
338 / 307 


3 94+0.27 
^■^ -0.27 

21+"-"' 
46.1+2] 
35.00!j5^ 
362 / 307 


4 91+U.26 
^■^^-0.20 

17+""*' 
"■^'-0^06 

25 3+" ' 
"■■'-0.9 

10 76+"" 

lU. /0„jj 

347 / 307 


4 89+"-2' 

^■°^-Q.28 

11+""*' 

-0^06 

15 1+"-^ 
6 03+" 26 

°-^"-0.26 
349 / 307 


4 61+"'" 
^■"'-0.43 

0.2 (fix) 

f. ^;+0.3 

2 66+""' 

^■™-0.16 
349 / 308 


9 5Q+1.00 
^■■'^-0.58 

0.2 (fix) 
2 0+"-2 
73+0 04 

'-'■'^-0.04 

333 / 308 


1 99+0.64 
'■^^-031 

0.2 (fix) 
3+" 2 

11 +""1 
"■"-0.01 

356/308 



*: Normalization of the apec component scaled with a factor SOURCE-RATIO-REG/fif from Table 3, 

Norm= souRCH-RAi'io-REG J- rt^„^JV7(4;r(l + z2)d2 ) x 10^20 ^^-5 arcmin"2_ where Da is the angular diameter distance to the source. 

t: Photon flux in units of 10~* photons cm"2 s^' arcmin"2 Energy band is 0.4 - 10.0 keV. 

Surface brighmess of the apec component scaled with a factor SOURCE-RATIO-REG and f2<, from Table 3. 



Table 6. The best fit parameters of \kT +PL model with statistical error 





18' 


-21' 


21' 


-24' 


24' 


-27' 




PL (r = 2.0) 


PL (r = 2.5) 


PL (r = 2.0) 


PL (r = 2.5) 


PL (r = 2.0) 


PL (r = 2.5) 


kT 
Z 
norm* 

r 


'■'■^^-(106 

14 S+^^s 
2.0 (fix) 


^ ^ ' -p.p6 

14 6+"^^ 
2.5 (fix) 


4.80!-' 

"■^ -0.3 
2.0 (fix) 


4 78+**-^'' 
27+"^i^ 
6 3 -^'^■^ 

V.J 

2.5 (fix) 


3 11+i^^f' 

0.2 (fix) 
1 Q+0.2 
^■^-0.3 

2.0 (fix) 


-'■^'■'-0.64 

0.2 (fix) 
1 9+0.3 

2.5 (fix) 


Flux UU 


1.02 


1.06 


0.99 


1.03 


0.99 


1.08 


xVd.o.f 


351/306 


351/306 


356/306 


355/306 


347/307 


348/307 



*: Same as Table 4. 



t: 10-40 keV flux of the power-law component, in 10 i* erg cm 2 s larcmin 2 



3.5. Constraint on Non-Thermal Emission 

We evaluated the flux of the possible non-thermal X-ray 
emission using three XIS spectra around the radio relic region 
(Fig. 3g, h, i). Kawano et al. (2009) gave the 3<t upper limit 
flux for a power-law component in the radio relic region to be 
1.1 X 10-^2 erg cm-2 s"' in 4-8 keV and 1.9 x 10"^^ erg cm ^ 
s~^ with extrapolation to 15-50 keV, respectively, assuming 
the photon index to be 2.0. We added a power-law compo- 
nent to the above ICM model as: apecuiB + wabs (apecuwH 
+ powerlawcxB + apecicu + power/awNon-themai)- The depen- 



dence on the photon index was examined by adopting 2.0 and 
2.5 in the fit. 

The results of the spectral fits are summarized in Table 6. 
We found no significant difference against the photon index, 
and aU three regions showed only upper limits for the power- 
law component. The derived upper-limit fluxes are f 10-40 kev < 
1.1 X 10"'^ erg cm"- s"' arcmin"^ at F = 2.5. To compare with 
the result by Kawano et al. (2009), we also derived the 4-8 keV 
flux as F4_8 keV < 5.2 X 10"^^ erg cm"-^ s~^ arcmin"-^ Using the 
solid angle covered by the radio relic (Qreiic = 122 arcmin^), the 
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Fig. 4. Radial profiles of (a) ICM temperature, (b) surface brightness, (c) deprojected electron density, and (d) gas pressure. Suzaku best fit values with 
statical errors are shown with black diamonds. Black dashed vertical lines show approximate radial boundaries of the west radio relic. Green and red 
dashed lines show typical systematic changes of the best-fit values due to change of OBF contaminants and NXB level. In temperature profile, gray 
dashed diamonds show the result of 2kT model fitting (see text). In surface brightness profile, blue and cyan dashed lines show the Galactic background 
and CXB emission, respectively. 



flux upper Umit in the 4-8 keVband is 5.3 x 10"'^ erg cm"^ s"\ 
which is consistent with the previous Suzaku results with XIS 
& HXD (Kawano et al. 2009). We note that the radial tempera- 
ture profile shows an increase just inside of the radio relic over 
the level of the central region by about 1 keV. The result of the 
spectral fit implies that this increase is not due to the presence 
of a hard power-law emission but rather showing the pure ICM 
heating caused by the shock. 

4. Discussion 

Suzaku performed 2 pointing observations of Abell 3376 
along its merger axis. The ICM temperature was found to be 
fairly flat at about 4 keV with a slow rise with radius in the 
region from the center to 0.6r2oo (1 Mpc). This is followed by 
a sharp drop in the radio relic region, which implies there is 
a shock front. We will attempt to evaluate the gas properties 



related with the shock (temperature, Mach number, magnetic 
field, and entropy) and discuss their implications. 

4.1. Temperature Profiles 

The ICM temperature is an important observational parame- 
ter which reflects the depth of the gravitational potential of the 
cluster. Their profiles also enable us to look into the history of 
the ICM heating and the growth of clusters. 

Previous Chandra, XMM-Newton results for other clusters 
indicated a significant temperature decline toward the outer re- 
gions. Recent studies with Suzaku showed temperature profiles 
of ICM to the virial radius (r2()o) for several clusters (George et 
al. 2008; Reiprich et al. 2009; Bautz et al. 2009; Hoshino et al. 
2010; Kawaharada et al. 2010; Simionescu et al. 2011). These 
profiles for relaxed clusters are approximated by a "universal" 
profile which declines by a factor of ~ 3 at r20() from the center. 
On the other hand, merging clusters show complex temperature 
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distribution depending on the mass ratio and collision geome- 
try. In the outer regions toward r2oo, radial temperature profiles 
show some excesses and sharp drops which are not seen in re- 
laxed systems (Akamatsu et al. 2011a; Akamatsu et al. 201 lb). 
However, very little systematic studies have been made about 
the ICM properties in the outer regions of merging clusters. 

We normalized the radial temperature profile by the average 
temperature determined within 0.4r2oo- Since A 3376 does not 
show a strong cool core, this should give a reasonable average 
temperature. The resulting scaled temperature profile is shown 
in Fig. 5 along with the previous Suzaku results (black crosses 
and diamonds: taken from Akamatsu et al. 2011a). Clearly, 
there is a marked difl'erence between the profiles of relaxed and 
merging clusters. The relaxed clusters (black crosses) gener- 
ally show a smooth decline from the center to the outer regions, 
consistent with the result of hydrodynamic simulations (Burns 
et al. 2010). Compared with these profiles, the A 3376 profile 
shows a large excess within 0.7r2oo before the sharp drop at 
0.7-0.8r2oo. 

Numerical simulation of merging clusters by Paul et al. 
(201 l)shows a shock feature very similar to the one in A 3376. 
This simulation after about 3 Gyr after a major merger shows 
an arc-like high temperature region whose width is about 200 
kpc, very similar to the observed temperature enhancement in 
Fig. 5. This suggests that A 3376 is still evolving after several 
Gyrs elapsed from the last major merger. 

The ACDM model predicts that clusters grow through sub- 
cluster mergers and matter accretion from large-scale struc- 
tures. Those dynamical events should give strong impacts 
on the ICM temperature structure, such as seen in A 3376. 
Since the ICM density is very low in the cluster outskirts 
(rie = 10"^ - 10"^ cm~^), the radiative cooling and conduc- 
tion times are expected to be very long with the cooling time 
largely exceeding the age of the Universe. We can naturally 
expect that the temperature structure carries important infor- 
mation about the history of the cluster growth. However, as 
mentioned above, there are no clear signs of dynamical evolu- 
tion in relaxed clusters. This indicates that the observed pecu- 
liar temperature structure in the merging clusters will be settled 
to the universal profile in a fairly short time scale, much shorter 
than the Hubble time. Since cooling time is too long, we have 
to consider other possibilities such as (i) ICM cooling by adi- 
abatic expansion, and (ii) ICM diffusion caused by pressure 
gradient. We will discuss these time scales in section 4.4. 

A fraction of the merger energy would be channeled into the 
acceleration of ultra-relativistic particles and amplification of 
magnetic fields. However, it is unknown that how these pro- 
cesses actually take place and distribute the merger energy to 
other forms in an efficient way. The turbulent motion can be 
detected through fine broadening by X-ray microcalorimeters 
such as with SXS instrument on ASTRO-H (Mitsuda et al. 
2010) and Athena, since X-ray calorimeters provide superior 
energy resolution by a factor of 20-30 better than that with the 
CCD instruments. The wealth of information will open a new 
window in our understanding of heating and particle accelera- 
tion caused by cluster merger shocks. 



4.2. Mach Number 

Recent studies with XMM-Newton and Suzaku showed 
sharp temperature drops at the region of the radio relic in 
A 3667 (Finoguenov et al. 2010; Akamatsu et al. 2011b). 
Finoguenov et al. (2010) reported a sharp edge in the sur- 
face brightness at the outer boundary of the NW radio relic in 
A 3667. Akamatsu et al. (201 lb) reported a significant jump in 
the ICM parameters (temperature, surface brightness, electron 
density, and pressure) across the same relic. The derived Mach 
number is Al ~ 2. 

We estimate the Mach number based on the Suzaku data in 
the same way as Akamatsu et al. (2011b). The Mach number 
can be obtained by applying the Rankine-Hugoniot jump con- 
dition, 

T2 _ + 14M^ - 3 1 

Ti ~ 16M2 ' C ~ 4M2 "^4' -* 

where C = ^ is the shock compression and subscripts 1 and 
2 denote pre-shock and post-shock values, respectively, assum- 
ing the ratio of specific heats as y = 5/3. Here, we assume 
the regions 21'-24' and 27'-31' to be the post and pre shock 
regions, respectively. Because the compression factor, C, ex- 
ceeds the value for the strong shock Umit, we only derived 
mach number from temperature jump. Table 7 shows the re- 
sultant Mach number, whose values are M - 2.94 + 0.77 based 
on the jumps of temperature. Compared with other clusters 
which also show shock fronts (A520, 1E0657-558, A2246, 
A 3667: Markevitch et al. 2005; Clowe et al. 2006; Russell et 
al. 2010; Akamatsu et al. 201 lb), the Mach number in A 3376 
is slightly larger than the other cases. 

We derive the compression factor, C - 7.8 + 4.7, which al- 
most exceeds the value for the strong shock limit, even though 
with a large error. The electron density estimated here is based 
on the assumption of spherical symmetry, which is not strictly 
correct when the shock propagates along the merger axis. Since 
the spherical assumption relates the observed flux to a shell- 
like volume with a larger line-of-sight depth than the case of 
a relic volume, the density in particular for the pre-shock re- 
gion is underestimated by a factor of roughly 1.5. In addition 
to this uncertainty, we checked the dependence on the centroid 
position of the spherical shock. We set the center of the sphere 
at the secondary X-ray peak which is closer to the radio relic. 
The resultant compression parameter is C = 5.01. This indi- 
cates that the determination of the centroid have a significant 
impact on the shock parameter, which is mainly resulted from 
the estimated line of sight depth of the shocked region. 

The shock front in A 3376 is located far from the cluster 
center (1.2 Mpc ~ 0.6 r2oo) just as in A 3667. The tempera- 
ture in the upstream region of the shock is ~ 1 keV, indicating 
the pre-shock sound speed to be ~ 520 km s Combining 
this with the shock compression C, we can evaluate the shock 
speed by fshock = C" • fss to be > 2080 km s"' assuming C = 4.0 
from the above discussion. The estimated shock speed well 
agrees with those in other clusters (1E0657-558: 4500 km s"', 
A520: 2300 km s"'), but twice higher than the A 3667 case 
(1360 ± 120 km s"'). The shock speed is large enough to ac- 
count for the radio reUc by a merger shock model (Takizawa & 
Naito 2000; Ricker & Sarazin 2001; Mathis et al. 2005). On 
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Table 7. Mach number estimation for A 3376 assuming 7 = 5/3. 



Parameter 


Pre shock 


Post shock 


Mach number 


kT |kcV| 


1.34 + 0.69 


4.68+0.48 


2.97+0.77 



the other hand, this ^shock seems too high considering the ICM 
temperature of 4 keV. We may need to consider a possibihty 
that the true gas temperature is much higher than the measured 
electron temperature. 

4.3. Entropy Profile 

Theoretical studies of the ICM (Tozzi & Norman 2001; Voit 
et al. 2003) predicts the entropy profile approximated by \ 
assuming that the gravitational energy of the accreting gas is ef- 
ficiently converted into the thermal energy. However, Chandra 
observations revealed that the entropy in the central regions of 
relaxed clusters often showed a flat profile (Cavagnolo et al. 
2009). The origin of this behavior is considered to be the feed- 
back of star formation and central AGN activities, which gen- 
erate extra entropy (Peterson & Fabian 2006). Entropy profile 
is thus sensitive to the non-gravitational energy injection and 
can be used to explore such a process even in the outer regions 
of clusters. 

We wiU define the "entropy" of the ICM hy K - kT ■ n^^'^. 
Fig. 5 shows the resultant entropy profile of A 3376. The slope 
shows a large deviation from the standard value of 1 . 1 and close 
to 0.4 in the radius range 0.'5-20.'0. The soUd black line shows 
K oc r^ '^, compared with the solid gray line indicating K oc r' '. 
Recent XMM-Newton study reported entropy profiles of 31 
clusters (Pratt et al. 2010), and the distribution of the slope 
indicated a peak value of 0.98. Therefore, the slope of 0.4 is 
very unusual and can be regarded as a transient value occurring 
in a certain phase of cluster mergers. 

The slope shows a steepening and flattening across the radio 
relic, which itself indicates the shock heating. A similar feature 
was also reported in A 3667 (Akamatsu et al. 2011b). The 
slope of A 3667 showed a marginal increase and a sudden drop 
across the radio relic region. In A 3376, the increase of the 
entropy slope is more significant, which suggests a stronger 
shock than in A 3667. This probably corresponds to the feature 
that the estimated Mach number in A 3376 is factor of 2 higher 
tiiatin A3667. 

4.4. Relaxation Time Scales 

As shown in the previous sections, A 3376 exhibits distinct 
ICM properties which are significantly different from those 
seen in the relaxed clusters. The density and temperature struc- 
tures around the shock region, along with the irregular mor- 
phology of this cluster, suggest that A 3376 is a young system 
regarding the relaxation of the ICM. The fact that most of the 
relaxed clusters now attain almost universal temperature profile 
indicates that the time needed for the ICM to reach thermal and 
dynamical equilibrium is much shorter than the cluster life. In 
this section we evaluate several relevant time scales based on 
the present observational quantities. 

We first estimate the time in which the shock front reaches 
the current position. Assuming that the shock has traveled the 
distance between the cluster center and the relic position with a 



constant velocity (v - 2000 km s~ ), the time required to prop- 
agate this length of ~ 1.5 Mpc is 0.32 Gyr. 

Next, we estimate the time scale to attain thermal equihb- 
rium after the shock heating. The shock heating should first 
act on ions rather than electrons, and then the thermal energy is 
transferred from ions to electrons through Coulomb collisions. 
The equilibration time scale between ion and electron is 3 or- 
ders of magnitude longer than the electron-electron timescale. 
Therefore, the cluster thermal time scale after the shock heat- 
ing is limited by the electron-ion time scale. According to 
Takizawa (1998), the electron-ion equiUbration timescale is 
given by 

where In A denotes the Coulomb logarithm (Spitzer 1956). 

In the central region, this gives 0.03 Gyr, which is about 10 
times shorter than the elapsed time after the shock passage. 
We calculated the elapsed time after the shock passage and the 
time for electron-ion equilibration fie as a function of radius. 
The resultant time scales are shown in Fig. 5 (c). Here, red di- 
amonds show fie, and gray dashed diamonds are the ionization 
time scale f, of Fe-Ka assuming net - 3x 10^^ cm"^ s. Solid 
line shows the elapsed time after the shock passage, assum- 
ing the shock propagation from iimer to outer regions. We can 
see that the elapsed time is shorter than fie outside of 0.6rshock 
corresponding to 17' (920 kpc). This suggests that the electron 
temperature is hkely to be lower than the ion temperature in the 
outer region including the radio reUc (Akahori & Yoshikawa 
2008; Akahori & Yoshikawa 2011). Those sign already re- 
ported another merging cluster RX J1347.5-1145 (Ota et al. 
2008). In the region that f, exceed the sound crossing time tsc 
( ti > tsc), there will be non-equilibrium ionization state. The 
future high resolution X-ray spectroscopy can reveal those dy- 
namical ionization states. 

FinaUy, we estimate the time scale in which the sharp change 
of the temperature disappears. Those high temperature re- 
gion will diffuse out due to the high pressure in the post- 
shock region and reach equihbrium in the sound crossing time 
(tsc = R/Vs), where R is the width of the high temperature re- 
gion and Us is the sound speed, respectively. Using the observed 
temperature (kT = 4.7 keV which corresponds to = 1 100 km 
s~^) and the width of the high temperature region {R = 320 
kpc), we can calculate the sound crossing time. The resultant 
time scale is tsc - 0.28 Gyr. 

In summary of this section, the observed features of temper- 
ature change, entropy slope and the time scale consideration all 
point to that the outer region (> 1 Mpc) of A 3376 is not reach- 
ing the equilibrium. In this sense, A 3376 is still a very young 
system in view of the cluster evolution. 

5. Summary 

We observed Abell 3376 with Suzaku XIS and derived radial 

profiles of temperature, surface brightness, electron density, 
and thermal pressure. The temperature and surface brightness 
shows remarkable jumps across the radio relic, located at ~ 1.5 
Mpc from the cluster center. We evaluated the Mach number 
using the Rankine-Hugoniot jump condition as = 3.0. The 
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Fig. 5. Radial profiles of normalized temperature, entropy and equilibration time scale, (a) Scaled projected temperature profiles compared with relaxed 
clusters (Akamatsu et al. 2011a). The profiles have been normahzed to the mean (within 0.4 raoo) temperature. The noo value is derived from Henry et 
al. (2009). Dotted line shows simulation result (Burns et al. 2010). Two gray dashed fines show standard deviation. Black crosses are for the relaxed 



clusters. Diamonds show merger clusters (black: A2142, blue: A3667, red: A3376). In entropy profile (b), black solid line shows K c 



given by Voit 



et al. 2003. (c) Ion-electron equilibration time scale. Horizontal scale is normahzed by the radius of the radio relic. Solid curve shows the time after a 
shock heating assuming a constant shock speed v = 2000 km s"' . Red diamonds show the electron-ion equilibration time Gray dashed diamond show 
the ionization time for n,t - 



main results on A 3376 are summarized as follows; 

• The ICM temperature is fairly flat at about 4 keV in the 
radius range from the center to 0.6r2oo ~ 1 ■ 1 Mpc, fol- 
lowed by a small enhancement near the radio relic to 5 
keV. 

• Across the radio relic region, the ICM parameters (tem- 
perature, electron density, and pressure) show significant 
jumps. The temperature drops from 4.5 keV to 1.0 keV, 
and the density and pressure both show an order of mag- 
nitude drop. 

• The estimated Mach number of the shock is Al = 2.97 ± 
0.77 based on the temperature jump, leading to the shock 
speed of yshock > 2000 km s"' . 

• The temperature structure observed in A 3376 is 
markedly different from those in relaxed clusters, sug- 
gesting that this cluster is a young system with the ICM 
properties strongly influenced by the merger shock. 

• The entropy profile is significantly flatter than the stan- 
dard prediction of r' ' and is approximated by r^*, ac- 
companied by slope changes across the radio relic. This 
indicates that the heating process is still going on in this 
region. 

• Based on the measurement temperature and electron 
density, the ion-electron relaxation time is longer than 
the elapsed time after the shock passage near the shock 
region. 

These results show that A 3376 is a dramatic merging clus- 
ter, with their outskirts still under a non-equilibrium condition. 
Because of the existence of the large scale shock front, future 
X-ray and SZ observations will provide new insight on the dy- 
namical evolution of clusters. 

The authors thank all the Suzaku team members for their 
support of the Suzaku project. We also thank an anonymous 
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(JSPS) Fellows (22-1582). 

Appendix 1. Stray Light 

To assess the effect of stray light from cluster core to out- 
skirts, we calculated photon contribution from each region to 
"On-Source" region using xissim (Ishisaki et al. 2007). We 
use y6-model (J3 - 0.40, = 2.'03) surface brightness profiles 
for input image and observed spectral information for input 
spectrum. Table 8 shows the resultant photon fraction detected 
in each region. The "On-Source" fraction shows the highest 
value followed by those in the adjacent regions, which is nat- 
urally understood considering the point-spread function of the 
Suzaku XRT (Serlemitsos et al. 2007). 

Appendix 2. ICM emission in the outermost region 

Since the validity of our detection of the ICM emission in 
the outermost region, r = 27' - 31', is important in this paper, 
we performed two additional checks. One is to look into the 
Suzaku spectrum by employing different background data, and 
the other is the analysis of the ROSAT pointing data for this 
cluster. 

To examine the adequacy of our background estimation 
process, we incorporated another background data which 
were taken at 4 degrees offset from A3376 on TX COL 
(ObslD: 404031010). We evaluated the background compo- 
nents (MWH and LHB) in the same fashion as described in 
Sec 3.2. The resultant best-fit spectral parameters are shown 
in Table. 9. Compared with the previous background estima- 
tion (Table 2), the temperatures and normalizations are con- 
sistent within the statistical errors, indicating that our back- 
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Table 8. Relative count contributions (%) due to PSF broadening of the Suzaku XIS mirror. Most counts come from the "On-Source" region. 
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Detector/ Sky 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(1) 0.'0-2.'0 


71.4 


26.2 


1.9 


0.4 


0.1 


0.0 


0.0 


0.0 


0.0 


0.0 


(z) Z.U-4.U 


Ll.i 


62.0 




1 'J 
1.3 


O.z 


0.0 


0.0 


0.0 


0.0 


0.0 


(3) 4.'0-6f0 


3.7 


25.1 


58.3 


12.2 


0.6 


0.0 


0.0 


0.0 


0.0 


0.0 


(4) 6.'0-9.'0 


1.2 


3.2 


18.7 


70.2 


6.4 


0.1 


0.1 


0.1 


0.0 


0.0 


(5) 12.'0-15.'0 


0.2 


0.3 


0.9 


10.6 


76.3 


10.1 


1.0 


0.3 


0.0 


0.2 


(6) 15.'0-18.'0 


0.1 


0.2 


0.0 


0.3 


11.7 


70.3 


15.9 


1.1 


0.2 


0.1 


(7) 18.'0-21.'0 


0.0 


0.2 


0.2 


0.2 


1.8 


13.4 


68.4 


14.7 


0.9 


0.2 


(8)21.'0-24.'0 


0.2 


0.2 


0.2 


0.1 


0.7 


1.3 


17.0 


65.9 


13.7 


0.6 


(9) 24.'0-27.'0 


0.2 


0.4 


0.3 


0.4 


0.5 


0.7 


1.6 


20.0 


66.3 


9.7 


(10) 27.'0-31.'0 


0.2 


0.3 


0.2 


0.4 


0.7 


0.3 


1.0 


2.5 


21.3 


73.1 



ground estimation is appropriate. Using the above background 
components, we analyze the outermost region (27' -31' annu- 
lus). The resultant best-fit values are kT - X.Tn^^^^^^ keV and 
norm=0.3 ± 0.2, which are consistent with the original one. 

Next, we derived the surface brightness of A3376 west di- 
rection using the ROSAT archival data (seqID: RP800154N00). 
We set the cluster center at (20*'12"3P,-56°49'12") and ex- 
tracted a region with a width of 6' to the west from the A3376 
center. We confirmed that the surface brightness dropped by 
an order of magnitude from r = 20' to 30' (Fig. 6). The bright- 
ness boundaries defined by the systematic error, mostly due 
to the background fluctuation, are indicated with blue his- 
tograms. Note that the ICM shows significant X-ray emis- 
sion around r = 30' even considering the systematic error 
(rp800154n00Jbkl.fits). These results indicate the validity of 
the Suzaku detection of the ICM emission in the outermost re- 
gion. 
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Table 9. Best-fit background parameters using TX COL 



[Vol. , 



fcr(keV) 


norm* (xlO 


fcr(keV) 


norm* (xlO ^) 




NOMINAL 0.08 (fix) 


^■-'-4.7 






296 / 247 



*:NormaUzation of the apec component scaled with a factor l/400;r. 

NoiTn=^j5jj J ngnndV/{4n{l + z^)D^) X lO"'* cm"^ arcmin"-, where Da is the angular diameter distance to the source. 




Fig. 6. Radial profile of the surface brightness. Red line show the ICM signal from A3376. Black line shows the background component. Two blue lines 
show the systematic en'ors considering fluctuations of the background component (±10%). 
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